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ABSTRACT 

We analyse INTEGRAL (between 2005 October and 2007 November) and 
RXTE (between 2007 June and 2011 March) observations of the accretion 
powered pulsar 4U 1907+09. From INTEGRAL IBIS-ISGRI and RXTE-PCA 
observations, we update pulse period history of the source. We construct power 
spectrum density of pulse frequencies and find that fluctuations in the pulse 
frequency derivatives are consistent with the random walk model with a noise 
strength of 1.27 x 10 -21 Hz s~ 2 . From the X-ray spectral analysis of RXTE- 
PCA observations, we find that Hydrogen column density is variable over the 
binary orbit, tending to increase just after the periastron passage. We also 
show that the X-ray spectrum gets hardened with decreasing X-ray flux. We 
discuss pulse-to-pulse variability of the source near dipping ingress and egress. 
We find that the source more likely undergoes in dipping states after apastron 
until periastron when the accretion from clumpy wind might dominate so that 
occasional transitions to temporary propeller state might occur. 

Key words: accretion, accretion discs - stars: neutron - pulsars: individual: 
4U 1907+09 - X-rays: binaries. 



1 INTRODUCTION 

The X-ray source 4U 1907+09 was discovered in early 1970s by the Uhuru survey (Giaconni 
et al. 1971). The system is an High mass X-ray binary (HMXB) that contains an X-ray pulsar 
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accreting material from its blue supergiant companion star. The pulsar has an eccentric (e 
~ 0.28) orbit around its companion and the orbital period of the binary system is ~8.3753 
days (In't Zand et al. 1998). 

The orbital profile of 4U 1907+09 exhibits two flares per orbit separated by ~ 0.5 orbital 
phase; a large primary at the periastron and the small secondary at the apastron (Marshall 
& Ricketts 1980, In't Zand et al. 1998). The presence of the two phase-locked flares had 
led to the suggestion that the compact object passes through a circumstellar disk of matter 
around the equatorial plane of a Be type stellar companion (lye 1986, Cook & Page 1987). 
However the position of 4U 1907+09 in the Corbet diagram (Corbet 1984) indicates that the 
companion is probably an OB type supergiant. Latest optical (Cox et al. 2005) and infrared 
(Nespoli et al. 2008) observations actually showed that the companion could be classifed as 
an 08 - 09 la supergiant with a mass loss rate of M = 7 x 10~ 6 M Q yr^ 1 , and a lower limit 
of distance ~5 kpc. 

The spin period of 4U 1907+09 was first measured as ~437.5 s using Tenma observations 
(Makishima et al. 1984). The pulse profile is known to be double peaked with variable 
amplitude. The profile is insensitive to the energies below 20 keV, whereas dramatic changes 
are present above 20 keV (In't Zand et al. 1998). Historical period measurements confirm 
that the pulsar had been steadily spinning down since its discovery until 1998 with an average 
rate of z> = -3.54 x 10~ 14 Hz s" 1 (Cook & Page 1987, In't Zand et al. 1998, Baykal et al. 
2001, Mukerjee et al. 2001). Afterwards, RXTE observations in 2001 showed that the spin 
rate was lowered by a factor of ~0.60 (Baykal et al. 2006). From INTEGRAL observations, 
it was reported that the spin period had reached to a maximum of ~441.3 s then, a torque 
reversal occurred and the source began to spin up with a rate of 2.58 x 10" 14 Hz s -1 after 
2004 May (Fritz et al. 2006). Recent measurements with RXTE presented in Inam et al. 
(2009a) and this paper have revealed that 4U 1907+09 has returned spin down trend with 
a rate of —3.59 x 10 -14 Hz s _1 , which is close to the previous steady spin down rate. This 
implies that another torque reversal should have been taken place before 2007 June. 

4U 1907+09 was identified as a variable X-ray source showing irregular flaring and dip- 
ping activities, ~20% observations of which were reported to be in dip state with no de- 
tectable pulsed emission (In't Zand et al. 1997). The typical duration of the dips was found 
to vary between few minutes to 1.5 hours. Variations in the X-ray flux were accepted to be 
the evidences of instability in the mass accretion rate where the dipping states are associated 
with the cessation of the accretion from an inhomogeneous wind of the companion star. 



The X-ray Pulsar 4 U 1907+09 3 

X-ray spectra of 4U 1907+09 were basically described by a power law with photon index 
~1.2 and an exponential cutoff at ~13 keV (Schwartz et al. 1980, Marshall & Ricketts 1980, 
Makishima et al. 1984, Cook & Page 1987, Chitnis et al. 1993, Roberts et al. 2001, Coburn 
et al. 2002, Baykal et al. 2006, Fritz et al. 2006). The continuum was found to be modified 
by highly variable Hydrogen column density (n H ) over the binary orbit, between 1 x 10 22 
cm -2 and 9 x 10 22 cm~ 2 as a consequence of inhomogeneous accretion via dense stellar wind 
(In't Zand et al. 1997). A narrow spectral line around 6.4 keV corresponding to Fe Ka 
emission produced by the fluorescence of matter surrounding the pulsar was also observed 
in the spectra. The detailed determination of Fe K emission complex has been recently 
reported using Suzaku observations, wherein Fe K(3 emission has been detected for the first 
time (Rivers et al. 2010). Observations with Ginga (Mihara 1995, Makishima et al. 1999) 
and BeppoSAX (Cusumano et al. 1998) exhibited cyclotron resonant scattering features 
(CRSFs) at higher energies. Line energies of the fundamental and the second harmonic 
cyclotron lines were found to be ~19 keV and ~39 keV respectively, implying a surface 
magnetic field strength of 2.1 x 10 12 G (Cusumano et al. 1998). 

In this paper, we present timing and spectral analysis of RXTE monitoring observations 
of 4U 1907+09 between 2007 June and 2011 March. Selected INTEGRAL pointing obser- 
vations are analyzed to cover the gap in spin period history. We describe the observations 
in Section 2. Pulse timing analysis is reported in Section 3. In Section 4, spectral results are 
discussed. In Section 5, we focus on the dipping states of the source. Finally, we summarize 
our results in Section 6. 

2 OBSERVATIONS 
2.1 RXTE 

One of the main instruments on-board RXTE is the Proportional Counter Array (PCA) 
which consists of five co-aligned identical proportional counter units (PCUs) pointed to the 
same location in the sky (Jahoda et al. 1996). The field of view (FOV) of PCA at full width at 
half maximum (FWHM) is about 1° and the effective area of each detector is approximately 
1300 cm 2 . PCA operates in the energy range 2-60 keV with an energy resolution of 18% 
at 6 keV. 

In this paper, we present timing and spectral analysis of 98 pointed RXTE-PCA ob- 
servations of 4U 1907+09 between 2007 June and 2011 March each with an exposure of 
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Table 1. Log of RXTE observations of 4U 1907+09. 



RXTE 


Time 


Number of 


Exposure 


proposal ID 


(MJD) 


observations 


(ks) 


93036 


54280 - 54825 


39 


74_ 1 <j,i; 


94036 


54839 - 55184 


21 


45. 3 b 


95350 


55208 - 55555 


26 


49.4 


96366 


55571 - 55638 


9 


47.2 



a Results from the first 30 observations were published before by Inam et al. (2009a). 
b Preliminary timing analysis results were presented before by §ahiner et al. (2011). 

~2 ks (see Tabled]). Preliminary timing analysis were already performed and corresponding 
spin-rate measurements were presented before by Inam et al. (2009a) and §ahiner et al. 
(2011). In this paper, we extend this analysis by using observations with a longer time span 
and performing a more detailed timing analysis as well as presenting our spectral analysis 
results. 

Although number of active PCUs during the observations of 4U 1907+09 varies between 
one and three, data obtained from PCUs and 1 are not appropriate for spectral analysis due 
to increased background levels as a result of the loss of their propane layers. Consequently, 
the only data products taken into consideration during spectral analysis belong to PCU 2 in 
order to avoid probable problems due to calibration differences between the detectors. The 
loss of propane layers of PCU and PCU 1 do not affect high resolution timing, therefore 
no PCU selection is done for the timing analysis. 

The standard software tools of HEAS0FT v. 6. 10 are used for the analysis of PCA data. 
Filter files are produced to apply constraints on the data such that; the times when elevation 
angle is less than 10°, offset from the source is greater than 0.02° and electron contamination 
of PCU2 is greater than 0.1 are excluded. Data modes examined for the spectral and light 
curve extraction are 'Standard2f ' and 'GoodXenon' modes respectively. Background spectra 
and light curves are generated by the latest PCA background estimator models supplied by 
the RXTE Guest Observer Facility (GOF) , Epoch 5C. 

2.2 INTEGRAL 

The INTEGRAL observations analysed in this paper are obtained from the INTEGRAL 
Science Data Centre (ISDC) archive. All publicly available pointing observations subsequent 
to the previous study of 4U 1907+09 (Fritz et al. 2006) are selected considering the good 
IBIS-ISGRI times to be above 1 ks. Selected observations, in which the source is in the FOV 
with an off-set angle smaller than 5°, are held on between 2005 October and 2007 November. 
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Figure 1. A sample template pulse profile (left) and its power spectra (right) in terms of harmonic number obtained from 
the observation on MJD 55586. 



Table 2. Pulse timing solution of 4U 1907+09 between MJD 54280 and 55600. 



Parameter 

Epoch (t ) 

Pulse Period (P) 

Spin frequency (v) 

v 

v 

ii 

"v' 

RMS residual (pulse phase) 



Value 



MJD 54467.61(6) 
441.2088(2) s 
2.2665002(8) X lO" 3 Hz 
-3.672(1) x 10- 14 Hz s- 1 
-1.497(5) x 10- 21 Hz s-2 
1.063(9) x 10~ 28 Hz s" 3 
-2.44(3) x 10~ 36 Hz s~ 4 
0.32 



The data consist of a total of 611 science windows (SCWs) (each ~3 ks) within revolutions 
366 and 623. 

The data products of IBIS-ISGRI detector on-board INTEGRAL are reduced for the 
analysis. Imager on Board the INTEGRAL Satellite (IBIS) is a coded mask instrument, 
which has a fully coded FOV of 8°.3 x 8°.0 and 12' angular resolution (FWHM) (Ubertini et 
al. 2003). INTEGRAL Soft Gamma-Ray Imager (ISGRI) is the upper layer of the IBIS in- 
strument which operates in the energy range 15 keV - 1 MeV with an energy resolution of 8% 
at 60 keV (Lebrun et al. 2003). The data reduction is performed by the software 0SA v . 7 . 0. 
The standard pipeline processing comprises gain correction, good-time handling, dead-time 
derivation, background correction and energy reconstruction. Images in two energy bands 
(20 - 40 keV and 40 - 60 keV) are produced from IBIS-ISGRI data with the use of an input 
catalogue consisting strong sources in the FOV: Ser X-l, XTE J1855-026, 4U 1909+07, SS 
433, IGR J19140+0951, GRS 1915+105 and 4U 1907+09. The background maps provided 
by the ISGRI team are used for background correction. Light curves are created by the tool 
II_LIGHT which allows high resolution timing. 
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Figure 2. Cycle counts of pulse phase and the best-fitting model (solid line). The residuals after the removal of 5th order 
polynomial (see Equation are given in the bottom panel. 

3 TIMING ANALYSIS 

3.1 Pulse Frequency Measurements 

For timing analysis, we use 1 s binned background corrected RXTE-PCA light curves of 
the source. These background subtracted light curves are corrected to the barycenter of 
the solar system. Then we correct the light curves for the binary motion of 4U 1907+09 
using the binary orbital parameters (In't Zand et al. 1998). Since 4U 1907+09 has many 
dips, we eliminate these dips from the light curve. Pulse periods for 4U 1907+09 are found 
by folding the time series on statistically independent trial periods (Leahy et al. 1983). 
Template pulse profiles are constructed from these observations by folding the data on the 
period giving maximum % 2 . The pulse profiles consist of 20 phase bins and are represented 
by their Fourier harmonics (Deeter & Boynton 1985). We present a sample template pulse 
profile and its power spectra in terms of harmonic number in Figure [TJ 

The pulse arrival times are obtained from the cross-correlation between template and 
pulse profiles obtained in each ~2 ks observation. We have been able to connect all pulse 
arrival times of RXTE observations in phase over a 1300 day time span. In the phase 
connection procedure, in order to avoid cycle count ambiguity, we construct the pulse arrival 
times for a time span where the maximum phase shift is less then 1. This time scale for 4U 
1907+09 is around 160 days. We divide total time span into 10 time intervals intervals each 
around 160 days and construct the pulse arrival times with respect to the best period in 
that time interval. Then, we align the slopes of pulse arrival times in the overlapping time 
intervals and construct pulse arrival times as presented in the upper panel of Figure [2j 




The X-ray Pulsar 4 U 1907+09 7 
We fit the pulse arrival times to the fifth order polynomial, 

5 H lift 1 

5<p = 5<p o + 6u(t - 1 ) + - toT (!) 

n=2 

where 5(j) is the pulse phase offset deduced from the pulse timing analysis, t a is the mid-time; 
5(f> is the residual phase offset at t ; 5u is the correction to the pulse frequency at time t ; 
4^ for n=2,3,4,5 are the first, second, third and fourth order derivatives of the pulse phase. 

The pulse arrival times (pulse cycles) and the residuals of the fit after the removal of the 
fifth order polynomial trend are presented in Figure El Table |2] presents the timing solution 
of 4U 1907+09 between MJD 54280 and 55600. It should be noted that we also obtained 
pulse residuals by performing pulse timing analysis using the timing solution parameters for 
pulse frequency and its time derivatives given in Table |2j 

In order to obtain pulse frequencies, we fit a linear model to each successive pairs of 
arrival times. Slopes of these linear fits lead us to estimate the pulse frequency values at 
mid-time of corresponding observations. In Figure [3] and Table we present the pulse 
period history of the source. In order to check the reliability of pulse arrival times, we also 
obtained the pulse periods using the slopes of pulse arrivals times obtained from 160-day 
time intervals. We found exactly the same pulse periods as obtained from the arrival times 
of 1300 days time span. Using the pulse periods, we obtained the pulse frequency value of 
i> = —3.672(1) x 10~ 14 Hz s _1 for a time span of ~1300 days . This is also consistent with 
the timing solution in Table |2j 

For the timing analysis of INTEGRAL observations, we extract 10 s binned background 
corrected IBIS-ISGRI light curves in the energy range 20 - 40 keV. We sample light curves 
from ~ 7 — 10 days time span and find the best pulse frequency by folding the light curve 
on statistically independent pulse frequencies. Then constructing master and sample pulse 
profiles, we obtain pulse arrival times as described above. From the slopes of pulse arrival 
times (5<f) = Su(t — t )) we attain the correction for pulse frequencies. We present pulse 
frequencies measured from RXTE-PCA and INTEGRAL IBIS-ISGRI observations in Figure 
[3] and Table [31 

3.2 Torque Noise Strength 

In order to see the statistical trend of the pulse frequency derivatives, we construct power 
spectrum of the pulse frequency derivatives. We use the Deeter polynomial estimator method 
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Figure 3. Pulse period history of 4U 1907+09. The measurements with RXTE and INTEGRAL in this work and Inam et al. 
(2009a) lie in regions A and B respectively and are separately plotted in the inset (see Table [3J. Measurements from Fritz et 
al. (2006) and Baykal et al. (2001, 2006) are in regions C and D respectively. Oldest measurements lie in region E (In't Zand 
et al. 1998, Mihara 1995, Cook & Page 1987, Makishima et al. 1984). 



(Deeter 1984) to the pulse frequency measurements. This technique uses the polynomial 
estimators instead of sinusoidal estimate for each time scale T. 

The power density spectra can be expressed by applying sampling function to the pulse 
arrival times (Deeter 1984; for the applications see Bildsten et al. (1997) and Baykal et al. 
(2007)). The power density estimator Po(f) is defined as J °° Po{f)df =< [y— < u >) 2 > 
where < v > is the mean pulse frequency derivative for a given analysis frequency. In 
order to estimate the power density, we first divide pulse arrival times into time spans 
of duration T and fit cubic polynomial in time. The observed time series is simulated by 
Monte Carlo techniques for a unit white noise strength defined as P„{f) = 1 and fitted to a 
cubic polynomial in time. Then the square of the third order polynomial term was divided 
into the value from Monte Carlo simulations (Deeter 1984, Cordes 1980). The logarithmic 
average of these estimators over the same time intervals is the power density estimate. 
This procedure is repeated for different durations T to obtain a power spectrum. We also 
calculated power density with quartic polynomial estimator and found consistent power 
density spectra with the cubic polynomial estimator. The frequency response of each power 
density and measurement noise estimates are presented at / ~ 1/T. In Figure HI we present 
power of pulse frequency derivatives (Pu(f)) per Hertz as a function of analysis frequency 
/. The slope of the power spectrum between 1/1300 and 1/75 d _1 is flat. This suggests 
that fluctuations in the pulse frequency derivatives are white noise and the pulse frequency 
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Table 3. Pulse period measurements of 4U 1907+09. 
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a The pulse periods are measured from INTEGRAL observations. 
b The pulse periods were presented before by §ahiner et al. (2011). 



fluctuations are consistent with random walk model. The noise strength is found to be 
1.27 x KT 21 Hz s- 2 . 

Random walk model in pulse frequency or white noise model in pulse frequency derivative 
are appropiate models for wind accretors e.g. Vela X-l, 4U 1538-52 and GX 301-2 (Deeter 
1981, Deeter et al. 1989, Bildsten et al. 1997). They have flat power spectra with white 
noise strength in the range 1CT 20 — 1CP 18 Hz s~ 2 . However, Her X-l and 4U 1626-67 are 
disc accretors with low mass companions which have shown pulse frequency time series 
consistent with the random walk model (Bildsten et al. 1997). Their white noise strengths 
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Figure 4. Power density of spin frequency derivatives of 4U 1907+09. The power due to measurement noise is subtracted from 
the estimates and shown independently by the cross symbols. 

are in the range 1CT 21 to 1CT 18 Hz s~ 2 . In these systems red noise in pulse frequency can 
not be ruled out (Bildsten et al. 1997). For disc accretors, Cen X-3 have red noise in pulse 
frequency derivatives, noise strength varies from low to high frequencies as 1CT 16 , 10~ 18 Hz 
s~ 2 (Bildsten et al. 1997). Power law index of power spectra in this system is ~ —1. This 
implies that, at short time scales, disc accretion is dominated and noise is less, on the other 
hand at time scales longer than viscous time scales there are excessive noise. 



4 SPECTRAL ANALYSIS 

4.1 Orbital Phase Resolved Spectra with RXTE-PCA 

The RXTE observations of 4U 1907+09 were proposed with a time distribution such that 
they spread over a variety of orbital phases, therefore they are very suitable for the inves- 
tigation of the variation in the spectral parameters through the binary orbit. During the 
analysis of PCA spectra energy range is restricted to 3 - 25 keV since the count statistics is 
poor beyond this range. The spectral analysis is carried out by using the XSPEC v. 12.6.0 
software. A systematic error of 2% is applied to handle the uncertainties in the response 
matrices and in background modelling (Wilms et al. 1999). 

4U 1907+09 is near the Galactic plane and the supernova remnant W49B, therefore an 
extra process is needed for a correct estimation of the background spectra. The exceptional 
dipping states of the source provide a good basis for background estimation since the count 
rates are consistent with the Galactic ridge emission (Roberts et al. 2001, Baykal et al. 2006). 
4U 1907+09 is detected to be in the dipping state in about 55 ks of the 216 ks RXTE-PCA 
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Figure 5. Variations of spectral parameters through binary orbit. The data points arc repeated for a cycle for clarity. From 
top to bottom; RXTE-ASM folded orbital profile, unabsorbed flux at 3 - 25 keV, Hydrogen column density and photon index 
plotted over orbital phase presented with 10 bins respectively. All uncertainties are calculated at 90% confidence level. The 
vertical dashed lines in all panels correspond to time of periastron passage vertical dotted lines in all panels correspond to time 
of apastron passage within 1 a. 

observations. The overall dipping state spectrum is used as an additional background, the 
procedure was explained before by Inam et al. (2009a). 

A total of 71 non-dip spectra of 4U 1907+09 are analysed for the orbital variation of 
the spectral parameters. Results of first 19 spectra were also given in the previous paper 
(Inam et al. 2009a). The basic model which consists of a power law (powerlaw) with a 
high energy cutoff (highecut) and photoelectric absorption (wabs) is successful for only 10 
spectra, where as an additional model component for the CRSF around 19 keV is required 
for the rest. We use cyclabs model in XSPEC (Mihara et al. 1990 and Makishima et al. 
1990). The fundamental line energy is fixed at 18.9 keV due to the statistical insignificance 
of individual PCA spectra with short exposure time ( ~2 ks) and the mean values of the 
depth and the width of the line are found as 0.4 keV and 1.7 keV respectively. These values 
are used to model CRSF since they agree within 1 a uncertainty with the previous CRSF 
parameters (Mihara 1995, Cusumano et al. 1998 and Makishima et al. 1999). Addition of a 
model component for the weak Fe emission line at 6.4 keV did not improve our fits, since 
this feature had lost its strength after subtraction of the diffuse emission from the Galactic 
ridge. 

Orbital variations of unabsorbed flux at 3 - 25 keV, Hydrogen column density (n#) and 
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photon index are plotted in Figure El which is an updated version of Figure 5 in Inam et al. 
(2009a). The variability of n# over the binary orbit is evident in the third panel of Figure 
|5j The base value of n# is ~ 2 x 10 22 cm -2 . It increases up to a value of ~ 10.5 x 10 22 cm -2 
just after the periastron passage (indicated by the vertical dash lines) and it remains at high 
values until the apastron, where it reduces to its base value again. 9 of our individual n# 
measurements exceed previously reported maximal value (~ 9 x 10 22 cm -2 , In't Zand et al. 
1997). Similar orbital dependence of was also reported for the source before (Roberts 
et al. 2001). One can see that there is no significant orbital variation of photon index in 
the bottom panel of Figure [5], the mean value is ~1.1. Other parameters that are consistent 
with being constant through the orbit are high energy cutoff and exponential folding energy 
with the mean values 12.7 keV and 9.6 keV respectively. 

In some accretion powered pulsars (e.g. Vela X-l (Haberl & White 1990)) n# is highly 
variable over the orbital phases and can range up to ~10 23 - 10 24 cm -2 . This variation 
is a probe to the density distribution of the accreting matter. The observed increase in 
absorption at periastron passage or eclipse ingress may be due to a gas stream from the 
companion trailing behind the pulsar (Haberl et al. 1989). In the case of 4U 1907+09, n# 
increasing to its maximum value after the periastron passage implies that the location of the 
absorbing material is the dense stellar wind of the companion star (Roberts et al 2001). Leahy 
(2001) and Kostka & Leahy (2010) modeled the absorption profile according to theoretical 
wind models and they have proposed that the most probable mechanism for 4U 1907+09 
is accretion from a spherical wind with equatorially enhanced dense spiral stream of gas 
around the companion star. 

4.2 X-ray Flux Dependence of Spectral Parameters 

To look for the X-ray flux dependence of spectral parameters of the source, we use spectral 
parameters and X-ray flux values obtained from 71 spectra that are presented in Section 
14.11 We sorted observations according to the X-ray flux values and spectral parameters of 
the observations with similar X-ray flux values are averaged to obtain a spectral parameter 
set as a function of X-ray flux. It should be noted that range of each flux value over which 
the spectral parameters are averaged is represented as x-axis error bars of the data points in 
Figure |6j We find that no spectral parameters except power law index show any correlation 
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Figure 6. Variation of power law index with X-ray flux from RXTE-PCA observations. 

with the 3-25 keV unabsorbed X-ray flux. Whereas power law index is found to be correlated 
with the X-ray flux (see Figure ED- 

Correlation between power law index and X-ray flux is an indication of spectral hardening 
with decreasing X-ray flux. Similar correlation has been seen in the accretion powered pulsars 
V0332+53 (Klochkov et al. 2011), 4U 0115+63 (Klochkov et al. 2011), 4U 1626-67 (Jain 
et al. 2010) and Vela X-l (Fiirst et al. 2009). On the other hand, there are also some 
accretion powered pulsars that exhibit spectral softening with decreasing X-ray flux: Her 
X-l (Kolochkov et al. 2011), Swift J1626. 6-5156 (Igdem et al. 2011), XMMU J054134.7- 
682550 (inam et al. 2009b), 2S 1417-62 (inam et al. 2004) and SAX J2103.5+4545 (Baykal 
et al. 2007). Both correlation and anti-correlation of power law index with X-ray flux might 
be considered as results of mass accretion rate variations and/or inhomogeneities in the 
companion's wind. 

5 DIPPING STATES AND PULSE-TO-PULSE VARIABILITY 

Since 4U 1907+09 is known to be a variable X-ray source, we analysed light curve of each 
RXTE-PCA observation to investigate pulse-to-pulse variations. The source is known to have 
dipping episodes that show no pulsations (In't Zand et al. 1997); therefore we searched for 
irregularities of pulsations due to dipping ingress and egress. Although short exposure time 
(~2 ks) of RXTE observations allow us to see at most about five successive pulses, we were 
able to detect a variety of rapid pulse shape changes. 

In Figured an example of a pulsating and dipping light curve of 4U 1907+09 is given in 
panel (a). Generally, the pulse-to-pulse X-ray intensity decreases just before dipping ingress 
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Figure 7. Sample set of 22 s binned light curves of different RXTE-PCA observations display the variations of pulse shape 
from pulse-to-pulse and the irregularities due to dipping activity. Time values are converted to phases (or time/pulse period) for 
arbitarly observation epochs. Explanations to each panel are given in Section[5] Observation IDs of panels from top to bottom are 
a:94036-01-19-00, b:96366-02-01-00, c:94036-01-05-00, d:94036-01-23-00, e:94036-01-13-00, f:93036-01-20-00, g:94036-01-03-00, 
h:93036-01-34-00, i: 93036-0 1-38-00 and j:95350-01-06-00. 
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(see panel (b)) and it increases gradually after the dipping egress (see panel (c)). However 
this is not the case for every dipping episode. The light curve in panel (d) is an extreme 
example, in which an X-ray flare follows a dipping episode and after a smooth decrease for 
two pulses, the source undergoes another dip again. Although the pulse shape is conserved in 
these primary examples, we find that it can also be totally disrupted by the dip, e.g. in panel 

(e) the pulse shape transforms to a single broad peak after a one pulse long dip and in panel 

(f) the attached wide peaks of the faint pulse are converted into detached sharp peaks after 
the dip. The beginning of the light curve in panel (f ) is also an example of unsuccessful dip 
which may be observed occasionally (see also In't Zand et al. 1997). In panel (g), single peak 
pulses are observable before and after the 1.5 phase long dip. In panel (h), after two broad 
pulses a dipping episode starts on phase 2.6 and continues until the end of the observation. 
A very rare condition is the total disappearance of only one peak of the pulse in panel (i), 
this example is the shortest dip we observe (between phases 2.9 and 3.4). We observe one 
pulse dips more commonly e.g. in panel (j) the dip between phases 1.8 and 2.8. 

Dipping states observed in Vela X-l (Kreykenbohm et al. 2008) have no identifiable 
transition phase in spectra before or after the dip. However, a recent study of GX 301-2 
(G6gii§ et al. 2011) reported a spectral softening starting before the dip. The photon index 
values reaching maximum during the peculiar dip of GX 301-2 resumes its normal values 
immediately after the dip. Fiirst et al. (2011) studied a dip spectra of GX 301-2 and found 
that although the pulses cease almost completely, nn does not show any significant variation. 
In the case of 4U 1907+09, direct spectral study of the dips with RXTE is not possible due to 
the diffuse galactic emission background, since the source is below the detection threshold 
and the dip spectra are similar to the Galactic Ridge spectra. However, it appears from 
Figure M that the source may be observed in dipping state in every orbital phase except 
between the phases 0.7 - 0.8 . Comparing this result with the spectral results of non-dipping 
observations (see Figure |5]); the absence of dips matches with the orbital phases when n# 
reaches to its maximum value just after the periastron passage (see lower panel of Figure |8]). 
Moreover, the possibility of dipping is less in orbital phases from 0.8 to 1.1 when compared to 
phases from 0.1 to 0.6; which is in anti-correlation with the n# values. The source frequently 
undergoes in dipping state during the times with minimum n# and the possibility peaks 
between the phases 0.4 - 0.5. It is important to note that this is not a strong conclusion, 
since direct measurement of the spectral parameters of the dips is not possible for us. 

Recently, strongly structured wind of the optical companion is accused of being the 
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reason of the dips in HMXB light curves. The clumpy wind of the companion, being highly 
inhomogeneous, creates regions of different density. Recent models have shown that the 
density of the wind material can change up to several orders of magnitude (Runacres & 
Owocki 2005). Therefore, a sudden decrease of X-ray luminosity is caused by the decrease 
in mass accretion rate. As the accretion rate decreases the Alfven radius increases and 
accretion stops when the source enters to propeller regime. Hence, the disappearance of 
pulsations can be explained by the passage of the pulsar through low density regions of the 
wind. The limit to the luminosity that turn on the propeller mechanism can be derived from 
the condition where the magnetospheric radius is equal to the co-rotation radius (Illarionov 
& Sunyaev 1975) and it is 6 x 10 31 erg s _1 for the case of 4U 1907+09. Therefore, the 
dipping states are candidate episodes for propeller regime since the count rates are below 
detection threshold. Similar propeller transition is suggested as the cause of observed dips 
in the accretion powered pulsars GX 1+4 and GRO J1744-28 (Cui 1997). 

The dipping episodes of 4U 1907+09 are known to have irregular time coverage, with 
varying duration. We investigate the distribution of dipping states through the orbital phase. 
About 46 of the 98 observations include dipping states. The percentage of exposure time 
spend in dipping is ~28% of the total exposure. In order to obtain the orbital distribution 
of dipping states, we calculate the percentage of the dip exposure to total exposure corre- 
sponding to each orbital phase with 0.1 phase resolution. The variation of the dip exposure 
percentage on orbital phase is plotted together with folded ASM light curve at orbital period 
at Figure EJ It is clearly seen that the occurrence of dips sharply decreases during the peri- 
astron passage and after the passage. Around the periastron passage, ASM count rate and 
therefore X-ray flux increases, whereas number of dip states decreases. It should be kept in 
our mind that a transient QPO had been found during one of the flares which suggests the 
formation of a transient disc (Int Zand et al. 1997). This implies that the number of occur- 
rences of dip states is less at the periastron passage and afterwards since the disc accretion 
is dominant mechanism. After the apastron (phases from 0.1 to 0.5), number of dip states 
increases due to the clumpy nature of accretion and occurrence of transient propeller regime 
episodes increases. The latter naturally explains the variable nature of the spin down. The 
combination of disc and wind accretion at different phases of orbit may lead to random walk 
nature of pulse frequency time series. 
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Figure 8. ASM light curve folded at orbital period (upper panel) and percentage of dip exposure times to the total exposure 
through the binary orbit (lower panel). The data points are repeated for a cycle for clarity. The vertical dashed lines correspond 
to time of periastron passage within 1 cr. It is evident that probability of observing a dipping state is lower during periastron 
passage and afterwards (between orbital phases 0.6 and 0.9). 



6 SUMMARY 

In this paper, we present analysis of INTEGRAL IBIS-ISGRI (between 2005 October and 
2007 November) and RXTE-PCA (between 2007 June and 2011 March) observations of the 
accretion powered pulsar 4U 1907+09. 

Using INTEGRAL and RXTE observations, we report new pulse period measurements 
and obtain an updated version of the pulse period history of Inam et al. (2009a) which is 
shown in Figure El Using these measurements, we construct power spectrum density of the 
pulse frequencies (see Figure HJ). We find that fluctuations in the pulse frequency derivatives 
are white noise and the pulse frequency fluctuations are consistent with the random walk 
model. The noise strength is found to be 1.27 x 10~ 21 Hz s -2 . We infer that in short time 
scales, there may be transient disc formation around the neutron star which causes random 
walk in the pulse frequency while in the long term the spin down rate is steady. 

From the X-ray spectral analysis of RXTE-PCA observations, we find that only Hydro- 
gen column density (n#) is significantly variable over the binary orbit, tending to increase 
just after the periastron passage (see Figure E]). This might indicate that the location of 
the absorbing material is the dense stellar wind of the companion star. We also study flux 
dependence of spectral parameters and find that power law index show correlation with the 
3-25 keV unabsorbed X-ray flux indicating spectral hardening with decreasing X-ray flux 
(see Figure |6]). This spectral variation might be related to mass accretion rate variations 
and/or inhomogeneities in the companion's wind. 

We also look for the irregularities of pulsations near dipping ingress and egress (see 
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Figure [7]). We propose that the disappearance of pulsations can be explained by the passage 
of the pulsar through low density regions of the clumpy wind which may lead to a turn-on 
of a temporary propeller state so that the X-ray luminosity decreases below the detection 
threshold. We are able to analyse orbital dependence of the dip state occurrences (see Figure 
|8]). We find that the source more likely undergoes in dipping state after the apastron until 
the periastron where transitions to temporary propeller state might occur due to accretion 
from clumpy wind. 
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